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2 '. Abstract 
O 

^p^' We study the possible impact of the neutrino oscillation which could be in- 

duced by a tiny violation of equivalence principle (VEP) for neutrinos emitted 
r~| , from supernova driven by gravitational collapse. Due to the absence of any 

significant indication of neutrino oscillation in the SN1987A data, we ob- 
' tain sever bounds on relevant VEP parameters ^ O(10~^^) for massless 

' or degenerated neutrinos and ^ O(10~^^) x [Am^/10~^ eV^] for massive 

^ neutrinos. 



I. INTRODUCTION 

Principle of equivalence, one of the fundamental bases in Einstein's general relativity, 
has been tested for many years with great accuracy. One can state the possible violation of 
equivalence principle (VEP) in several ways. For massive objects, existence of VEP implies 
that gravitational mass is not equal to inertial mass. For macroscopic objects, by the 
torsion balance experiment, Eotvos and his collaborators [|l| obtained the stringent bound 
?7 < 3 X 10^^, where 1] is the parameter which measures the difference between gravitational 
and inertial mass. Forty years later, Dicke et al. P| improved the bound to r/ < 3 x 10^^^ by 
removing some systematic uncertainties making use of the influence of Sun's gravitational 
field which could produce a torque with 24 hour periodicity in the experiment in the presence 
of VEP. With similar apparatus, Braginski and Panov 0] obtained the bound r] ^ 0.9 x 10~^^. 
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Equivalence principle has been tested also for microscopic objects. Using the neutron 
free fall refractometry, the bound 77 < 3 x 10~^ is obtained |^ for inertial and gravitational 
mass difference of neutrons. 

For particles like photon and neutrinos, one can not perform such a comparison of inertial 
and gravitational masses. One can use, instead, another parameter, ^7, in order to test 
VEP, which measures the difference of the gravitational coupling between the two particles 
involved, d-^pa = {Gp — Ga)/{G/^ + Go)- A mild bound \5'^'yv\ ^ 10~^ for violation of 
the equivalence principle for photons and neutrinos have been derived by using the 
small arrival time differences between photons and neutrinos from the supernova SN1987A 
observed by Kamiokande and 1MB detectors 0. 

For neutrinos, more sensitive test of VEP can be performed by using neutrino oscillation 
which can be induced by VEP even if neutrinos are massless Based on this proposal. 



many works have been performed ||9|-p3|. In Ref. ||2J] such VEP induced neutrino oscillation 
was shown to be phenomenologically equivalent to neutrino oscillations induced by a possible 
violation of Lorentz invariance 



Some theoretical insight on the type of gravitational potential that could violate the 



weak equivalence principle can be found in Ref. [26|. Some discussions on the possibility 



that violation of equivalence principle comes from string theory are found in Ref. I^^- A 
discussion on the departure from exact Lorentz invariance in the standard model Lagrangian 
is developed in Ref. [p8[] . 

A possibility to test VEP by the future long-base line neutrino oscillation experiment 
has been discussed in Ref. |19[ whereas some bounds from the current accelerator based 
neutrino oscillation experiments were obtained in Ref. |^. Recently such neutrino oscilla- 
tion mechanisms have been investigated [ll7| , p!8| in the light of the experimental results from 
Super-Kamiokande (SK) on the atmospheric neutrino anomaly, obtaining stringent limits 
for the — > Vr channel. 

The possibility of solving the solar neutrino problem (SNP) by these gravitationally in- 
duced neutrino oscillations have been investigated in Refs. [§-0. On the other hand, pos- 
sibilities to constrain VEP parameters by using astrophysical neutrino sources are discussed 
in Ref. [|1^] for solar neutrinos and in Ref. for very high energy neutrinos from AGN. 
The possible effect of neutrino oscillation induced by the violation of equivalence principle 
(VEP) as an explanation of the pulsar kick velocity has been considered in Refs. p^ , |23[| . In 
Ref. |Tl[], it was discussed that supernova (SN) explosion dynamics as well as heavy element 
nucleosynthesis can be significantly affected by VEP induced neutrino oscillation. 

In this paper, extending the discussions in Ref. |TT[], we consider the possible impact of 
VEP induced neutrino oscillation for SN neutrinos, in particular for the observation of 9e 
signal. In this work, we do not consider the possible effect of neutrino oscillation induced by 
VEP on heavy element nucleosynthesis in SN |]2D| as the bounds on VEP parameters from 
this effect is expected to be weaker. We show that one can try to constrain the relevant 
VEP parameters using the observed anti-electron neutrinos events from supernova SN1987A 



following the same argument used in Refs. p0|j31[]. In Sec. II, we discuss the formalism we 



use in this work. In Sec. Ill we describe some features of SN neutrinos and discuss how to 
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use supernova neutrino data to constrain oscillation parameters. In Sec. IV we present our 
results. In Sec. V we give conclusions. 



II. NEUTRINO PROPAGATION IN THE PRESENCE OF VIOLATION OF 

EQUIVALENCE PRINCIPLE 



Here we describe the framework we will use in this work. Let us assume that neutrinos 
of different species are not universally coupled to gravity, or neutrinos of different flavor 
propagate with different gravitational red shift in a given gravitational field. In such case, the 
weak interactions eigenstates l^w) and gravitational eigenstates |z/g), the bases with which 
gravitational part of the energy is diagonalized, can, in general, be different, l^w) 7^ \^g)- 
If neutrinos are massive, the mass eigenstates \i'm)) can be different from both of these two 



states, \h'M) 7^ \i^w) 7^ \^g)- For two neutrino system, Ve 
such a general situation as follows. 



/i or r), we can describe 
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where 9g is the mixing angle due to the presence of VEP which is in general different from the 
usual flavor mixing 9m between weak and mass eigenstates, 9g 7^ 9m, implying {um) 7^ Wg)- 
For simplicity, let us consider, the system of two neutrino flavors. The evolution equation 
of massive neutrino for one generation without flavor mixing in the gravitational field in the 
supernova can be described by. 



.d_ 
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(2) 



where the gravitational potential 0(r) describes the red shift of neutrinos and it is given, up 
to the first order in Newton's gravitational constant G, assuming the spherically symmetric 
metric and taking the center of the star as origin of the coordinate, as follows. 



(r) = -G 



'M(r) 
— — H 




r 





+ 



where 



M(r) 



dr'Anr' p{r' 



(3) 



(4) 



is the total mass contained in the volume with radius r, and (psc is the background potential 
which can be taken as the one coming from the local Super-cluster [B^. For definiteness, in 
this work, we take (psc = —3 x 10~^, the value estimated in Ref. |3^ . 

It is straight forward to generalize the evolution equation for more than one generation 
with mixing. Since, as we will discuss in detail in the following sections, we are most 
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interested in the oscillation effect for anti-electron neutrinos, Ue, from now on we will consider 
only Ue — i^x{x = ^i or r) system unless otherwise stated. The evolution equation for massive 
anti- neutrinos in the presence of VEP can be written as follows 



.d 
'di 







' -Ve - Vg cos 2eG - A cos 26*4^ 


|(Vbsin2^G + Asin2^Af)" 






_ 




i(l^Gsin2^G + Asin2^M) 










(5) 



where Ve = ^J^CppYelvrij^ is the standard matter potential, A = Ara^ /2E, 6m is the usual 
flavor mixing induced by mass which can be different from 9g and Vg{t) = 26'jE(f){r). For 
neutrinos system, z/g — Ux, the same equation hold by replacing —Ve by Ve- 

We assume that as a reasonably good approximation the density profile in the relevant 
region of the star takes the following form. 



P{r) = Po 



r 



(6) 



where m = 5—7 in the inner part and m = 3 in the outer layer of the star. For definiteness, we 
take Po = pi = 8x 10^"^ g/ cm'^ and Rq = Ri = 10^ cm for the inner part and po = P2 = 4 x 10^ 
g/ cm^ and Rq = R2 = 10^ for the outer part of the star. The radius which defines inner/outer 
part has been taken to be the distance at which both definitions merge, Ru = 4.5 x 10^ cm. 
Under above assumptions, 0(r) can be explicitly computed as. 
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12, 



(7) 



+ (psc for Ru < r, (8) 

where Mi = AirpiRf/S and M2 = Mi(i?i/i?2)^- In Fig. 1 we show the profiles of 
Ve + A COS 26m and Vq in eV for some values of mass and mixing parameters, and VEP 
parameters, as a function of distance from the center of the star. 

A resonantly enhanced conversion can occur if the the resonance condition as well as the 
adiabaticity conditions are satisfied. The resonance condition is given by, 

Ve = -Vg cos 2eG - A cos 26 m, (9) 

Note that depending on the sign of Vg and A, resonance condition can be satisfied either 
for neutrino or anti-neutrino. 

The adiabaticity parameter /t(r) at position r can be defined as 

(10) 

where AH{r) is the energy splitting between two energy eigenvalues of the Hamiltonian, 
AH{r) = — Emi- Explicitly, K,{r) can be written as follows. 
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where io = ^-KniN / \/2G p pYe is the matter refraction length, ^l, = A-nE j /S.rn? is the oscillation 
length in vacuum without VEP and l^J^ is defined as, 



G 



^sin2gG + f sin 2(^0- ^m) 
sin 2610 + ^sin26'M 



(12) 



dlnV^/dr| ^ and £g = IdlnVc/drl ^ are the density and gravitational scale 

(13) 



where (.p 

height, respectively. Mixing angle in matter 9m in the presence of VEP is given by, 

2Hi2 A sin 26 m + Vg sin 26*0 



tan 29r, 



if 22 - Hii K + A COS 26 M + Vg cos 26 g 



For the case where resonant conversion occur, at resonance point, k takes somewhat 
simpler form. 



, (A sin 26m + Vcivres) sin 26g) 



^ + ^cos26g 



(14) 



from which it is possible to derive the known cases, when Am? — or ^7 ^ 0. 

We compute the survival probability of at the surface of the outer region of the 
supernova , by using the following formula. 



pSN ^ p 

ee 



SN t 



Pr) COS 2^,- COS 2^ 



(15) 



where 6i and 6f are the initial and final mixing angle which can be taken as cos2ftj = 1 and 
cos 26 f = cos26'o where 6q is defined as 



tan2^^n 



Asin20A/ + V^"" sm26G 
A cos 2^A/ + V'g''" cos2^g' 



(16) 



The crossing probability Pc is computed as [^ , 

exp(-f kF) - exp(-f kF/ sin^ ^0) 



1 -exp(-fKF/sin2 60) 



(17) 



where, as an approximation, following Ref. [33|, we use F = \1 — tan26'o| which is valid 
for exponentially decreasing density profile. When does not undergo resonance, as an 
approximation, we use the same formulas but with k estimated at the resonance point for 
channel. 



We compute the final survival probability at the detector as follows ^4 
Pnii^e) 



SN 



sm' 60 + Piil ~2P^ 

~r -If 



COS 26^0 



(18) 
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where P^^ is the averaged survival probabihty of Ue at the surface of supernova, computed by 
taking into account the energy spectrum (see eq. (p3|)) and is the conversion probabihty 
of z/2 — > after crossing the Earth. We define new eigenstates, vi and V2i which are neither 
the mass nor gravitational eigenstates in the presence of mass and VEP, by using the mixing 
angle defined in eq. (p!6D as follows, 



(19) 
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Pcoh is the "coherent" part of the probability which is given by, 
2v/P,r (1 - Per)Pe^(l - Pi) cosiasN + + 



p 



coh 



(20) 



where P^/^ and P^ denote the conversion probability of Ue — ^ i^i in the supernova and in 
the Earth, respectively, asN is the relative phase difference between the ui and z/2 states at 
the surface of the supernova, a^; is the phase difference between the amplitude of z/i — >• 
and z/2 Vf. after crossing the Earth, and 



a 



-A- 



; sin 4:6g sin 46 



M 



(21) 



where L is the distance between supernova and the Earth. In our analysis we will use L = 50 
kpc for the distance between SN1987A and the Earth. We note that Pcoh could be important 
if sum of these phases is not so large. If one of them is very large, then the coherent term will 
be averaged out, after we take into account the energy spread of supernova neutrinos. The 
largest contribution is coming from a and when a is small, the other phases are expected 
to be very small compared to a, and therefore, we can simply neglect other two phases asN 
and aE, as a good approximation, to compute Pcoh- 



Assuming that the density is constant, can be computed as follows ^0 

Ue) = sin^ 6*0 + sin 26'® sin 2(6'© — 6*0) sin^ > 



P 



2e 



P^(^^2 



where i is the distance traveled by neutrino inside the Earth, 6'® and 
and the neutrino oscillation length, respectively, in the Earth. 



(22) 

are the mixing angle 



III. POSSIBLE IMPACT OF OSCILLATION FOR SUPERNOVA NEUTRINOS 

Let us consider the influence of neutrino oscillation between z/g and z/^. (x = yU or r) for 
SN neutrinos. From the view point of observation, influence on z/g signal is most important 
because the cross section for absorption reaction UeP —>■ e~^n is much larger than elastic 
scattering processes z/^-e" z/^-e" (x = e,fi,T). In fact it is considered that observed events 
at Kamiokande-II and 1MB detectors from SN1987A are induced by z/g. 

The energy spectrum of supernova neutrinos can be described by "pinched" Fermi-Dirac 
distribution where pinched form can be parametrized by the "effective" chemical potential 



rji as follows [35 
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(23) 



where Tj and rji are the temperature and effective chemical potential, respectively, for i-th. 
neutrino species. Because of the different interaction rates of neutrinos in the proto-neutron 
star, it is expected that there exist the following hierarchy relation between temperatures, 

T.^<Tp^<T,^_, =.Tp^,,. (24) 

Due to such temperature difference, average neutrino energies in typical numerical SN sim- 



ulations, which are approximately related with temperatures as (E,^) ^ 3T,y, are p6|-[38[] 



{E,J ^11-12 MeV, 
{EpJ ^14-17 MeV, 
(^M.))^(^M.))^ 24- 27 MeV. (25) 

Such energy hierarchy is of crucial importance for our discussion as we will see below. 

If neutrino oscillation between z/g and {x = fJ',T) occurs in SN and/or in the space 
between SN and Earth, expected neutrino flux at the Earth can be given by, 

F,^ = il-p)Fl+pFl, (26) 

where p is the permutation factor, or the energy averaged conversion probability for ^ i/^ 
which can be computed by convoluting the probability in eq. ([T8|) in the previous section 
with the Fermi-Dirac energy distribution in eq. (|23D, and F^ denotes the time integrated 
original neutrino flux in the absence of oscillation. For simplicity, we will set the "effective 
chemical potential, rji, in eq. (^3]), equal to zero in this work. 

Because of the energy hierarchy in eq. (|25|) , large oscillation between z/g ^ z/^,r would 
make the observable energy of z/g larger than the theoretical predictions. However, the 
observed SN1987A neutrino data by Kamiokande and 1MB detectors imply rather smaller 
average energy, which are in the lower side of the theoretical expectations, indicating that 
there was no significant oscillation between z/g and z/^^^. Based on this argument, it was 
obtained in Ref. [0 the bound on permutation parameter as p < 0.35 at 99 % C. L., which 



disfavors large oscillation between Vf, and far- The same argument was used in Ref. |31| and 



most recently in Refs. PU|-^, to constrain large oscillation for neutrinos from SN1987A. 

IV. BOUNDS ON VEP PARAMETERS FROM SUPERNOVA 

In this section we discuss our results which are obtained using the formulas described in 
the previous section. 
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A. massless neutrinos case 



First let us discuss the case without neutrino mass, or the case where masses are de- 
generated. In this case, the relevant formulas can be obtained by putting Am^ — in the 
previous section. In order to cover all the physical case, we fix ^7 to be positive but consider 
9g in the range from to 7r/2, analogous to the case of usual MSW effect Here we try 



to constrain the VEP parameters from the observation of the supernova SN1987A neutrinos 



by Kamiokande and 1MB following the argument discussed in Refs. ||3yj31[ . 

In Fig. 1^ we show the iso-contours of permutation factor for anti electron neutrinos from 
supernova SN1987A on the tan^ Oq — ^7 plane. For parameters in the inner region of these 
curves, the neutrino conversion is efficient so that the energy spectrum of Ve gets harder 
which would lead to worse agreement of the data with the predictions from SN simulations. 



If we use the criteria adopted in Ref. [0, shaded region could be disfavored at 99 % C.L. 
In the same plot, for the purpose of comparison, we also show the region (inside the dashed 
curves) which can be excluded by the laboratory experiments obtained in Ref. |2l| . 

We note that for ^7 larger than 10~^°, strong VEP induced resonant conversion P], 
which is similar to the well known MSW effect, can occur inside the supernova, and the 
permutation parameter can be larger than 0.5 even for very small values of gravitational 
mixing angle Oq- We see that significantly large parameter region can be disfavored which 
can not be tested by laboratory experiments (see dashed curve). 

On the other hand, for ^7 smaller than ~ 10^^°, the dominant effect is the VEP induced 
vacuum like-oscillation with the probability, 

P(z?e ^ z/e) = 1 - sin^ 260 si\i^[5-i<pscEL]. (27) 

In this case the permutation parameter can be large if the mixing angle 9g is large and 
57 ^ 10~^^ and therefore, such parameter region could be constrained. The reason why we 



can obtain such stringent bound on ^7, far below than the existing experimental limit |^ 
57 ^ lO~^^[|0sc|/3 X 10~^], is simply because of the fact that the distance between the 
observed supernova and Earth is much larger (about 50 kpc, for the case of SN1987A) 
compared to the baseline of laboratory neutrino oscillation experiments or even compared 
to the baseline for solar neutrinos, Sun-Earth distance. 



B. massive neutrinos case 

Next let us discuss the case where neutrinos are massive. Let us consider how the result 
presented in the previous section can be modified due to mass and usual fiavor mixing. For 
the sake of demonstration, we take some particular sets of mixing parameter implied by the 
MSW solutions to the solar neutrino problem ||^. Basically, the main effect is that the 



sensitivity to VEP is weakened due to the presence of usual mass and mixing terms. It is 
easy to see this. When the magnitudes of mass and usual flavor mixing terms are lager than 
those by VEP, oscillation effect is essentially induced or dominated by the usual mass and 
flavor mixing, and vice versa, as they mainly govern the evolution of neutrinos. 
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In Fig. 3 we show our results obtained by assuming, particular set of the mass and 
flavor mixing implied by the large mixing angle (LMA) MSW solution to the solar neutrino 
problem, Am^ = 3.2 x 10"^ eV^ and tan^OM = 0.33 in addition the VEP effect. For 
^7 smaller than ~ O(10~^^) the effect of neutrino conversion is essentially dominated by 
the usual mass induced oscillation. We note that new disfavored region appears for ^7 ~ 
10^^^ — 10"^'' for any small values of 6g- This is due to the combined effects of the usual mass 
and VEP induced terms in the neutrino evolution Hamiltonian. In these parameter regions 
of 57 around ~ O(10~^^), resonant conversion is mainly satisfled by the VEP parameters 
whereas the adiabaticity condition is met due to the presence of mass induced mixing term 
or |Asin2^M, even if 9g is zero. 

In Fig. 4 we show the similar plot but for the small mixing angle (SMA) MSW solution 
with Am^ = 5.0 x 10~^ eV^ and tan^ 6m = 5.8 x 10~^. In this case, unlike in Fig. 3, no new 
disfavored region appears for ^7 ~ 0(10^^^). This is because for these values of ^7 with 
small ^ 1; the adiabaticity condition, which was satisfied in the case of LMA solution, 
is not satisfied anymore since |Asin26'Af is much smaller than the case of LMA. 

Some comments are in order. We showed in Fig. 3 and 4 some parameter region which 
can be disfavored by SN1987A data assuming the mass and mixing suggested by LMA and 
SMA solutions to the solar neutrino problem. We should note, however, that if ^7 is much 
larger than ~ 10~^^, solar neutrino fiux can be significantly affected by the VEP effect 



and hence LMA and SMA solutions must be affected accordingly, as discussed in Ref. [15 



Therefore, the parameter regions allowed in Fig. 3 and 4 by the SN1987A data are not 
necessarily implying the parameter region allowed by the solar neutrino data. 



V. CONCLUSIONS 

We have studied the impact of the neutrino oscillation which could be induced by a 
tiny breakdown of the equivalence principle. We have obtained regions of parameters which 
are disfavored by the SN1987A neutrino data and showed that supernova can prove some 
region of parameter space which can not be tested by laboratory experiments. We show 
that supernova neutrinos can be sensitive to disfavor very tiny values of ^7 as small as 
O(10~^^), which is far below the laboratory bound, if neutrino mass can be neglected. We 
also showed that much smaller gravitationally induced mixing angle, tan^ 6g ^ 10~^, which 
can not be tested by the laboratory experiments, can be proved by SN neutrinos. For 
massive neutrinos, we showed that the sensitivity to these bounds will become worse. The 
effect of VEP in oscillation is lost or weakened when mass induced terms become larger 
than VEP terms. Let us note that, compared to the bounds on VEP obtained by using the 
macroscopic objects, rj < 10~^^, mentioned in the introduction, these bounds we obtained 
are much more stringent. 

Let us finally try to compare the sensitivity to VEP parameter for neutrinos from su- 
pernova with that for neutrinos from other astrophysical sources. With solar neutrinos, it 
was discussed that effect of VEP can be tested for ^7 as small as ~ 10^^° for massless 
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neutrinos [T^ and for ^7 as small as ~ 10^ — 10~ for massive neutrinos if we assume 
the MSW solutions to the solar neutrino problem |15|. Much better sensitivity can be ob- 
tained if we can observe very high energy neutrino from AGN []TB|. Due to much higher 
energy and lager distance, it is discussed in Ref. [16| that one can test VEP for ^7 as small 
as 10"'^^ if neutrinos are massless or degenerated neutrinos and ^7 ~ 10^^^ x (Am^/leV^) 
for massive neutrinos. However, let us stress that in the absence of such ANG neutrino 
data, currently, SN1987A data seem to provide better bound on VEP parameters than any 
existing terrestrial experiments. 
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FIGURES 




r (cm) 



FIG. 1. Voir) are plotted as a function of the distance from the center of the star for various 
values of ^7 by dotted, dot-dashed, short-dashed, long-dashed curves. Also Fmot('') + A cos 20 is 
plotted for massless (or degenerate) case and E = 10 MeV, Am^ = 10~^ eV^ and sin^ 29m = 0.9. 
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FIG. 2. Iso-permutation factor contours for anti electron neutrinos from supernova on the 
tan^ — ^1 plane. The curves correspond, from inside to outside, p = 0.9,0.75,0.35,0.23,0.1. 
The shaded region corresponds to p > 0.35. Inside the dashed curve indicates the excluded region 



obtained from laboratory experiment obtained in Ref. |21]. 
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